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PREFACE 
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(a  subsidiary  of  Sverdrup  & Parcel  and  Associates,  Inc.),  contract 
operator  of  AEDC,  Arnold  Air  Force  Station,  Tennessee.  The  research 
was  done  under  ARO  Project  Nos.  V43A-07A  and  V43A-16A.  The  author  of 
this  report  was  M.  0.  Varner,  ARO,  Inc.  The  manuscript  (ARO  Control 
No.  ARO-VKF-TR-76-95)  was  submitted  for  publication  on  August  20,  1976. 

The  work  of  Mr.  D.  E.  Boylan  in  obtaining  the  experimental  data 
in  Tunnel  C and  Dr.  J.  C.  Adams  in  providing  the  Lubard  hypersonic 
viscous  shock  layer  (HVSL)  theoretical  results  is  acknowledged. 
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1.0  INTRODUCTION 

The  total  temperature  probe  currently  used  In  the  von  Karman  Gas 
Dynamics  Facility  (VKF)  at  AEDC  is  of  the  single-shielded  type  as 
developed  by  Bontrager  (Ref.  1).  He  considered  the  cumulative  effects 
of  radiation  loss  to  the  walls,  velocity,  and  conduction  in  defining 
the  optimum  probe  configuration  for  a given  flow  environment.  Probe 
calibration  was  accomplished  by  defining  the  relationship  of  actual- 
to-measured  total  temperature  in  the  free  stream  at  Mach  8 and  10 
with  the  ratio  of  estimated  entrance  length  of  underdeveloped  flow 
to  a characteristic  probe  length.  An  experimental  parametric  analysis 
revealed  an  optimum  probe  entrance  length-to- internal  diameter  ratio 
of  two  and  probe  entrance-to-vent  area  ratio  of  two  or  greater  over 
the  range  of  variables  investigated. 

Recent  use  of  Bontrager 's  method  for  total  temperature  probe 
correc"  has  pointed  out  a shortcoming  of  his  approach.  As  long 
as  the  calibration  and  corrections  are  made  over  approximately  the 
same  Mach  number  and  probe  Reynolds  number  range,  the  method  appears 
to  give  consistent  results.  Use  of  an  extrapolated  calibration 
curve  to  correct  total  temperature  measurements  in  an  uncalibrated 
Mach  number  and/or  probe  Reynolds  number  range  has  led  to  inconsistent 
results  when  compared  to  theoretical  predictions. 

Others,  for  example,  Winkler  (Ref.  2)  and  Voisinet,  Lee,  and 
Meier  (Ref.  3),  have  included  a Mach  number  effect  based  on  the  total 
temperature  recovery  of  cylinders  in  compressible  flow.  This  approach 
as  applied  to  a somewhat  different  probe  design  appears  to  give  an 
improved  correction  to  the  total  temperature  measurements.  Apart  from 
the  correction  for  Mach  number  effect,  the  inaccuracy  still  exists  in 
extrapolating  a calibration  curve  based  on  tunnel  free-stream  conditions 
to  points  outside  of  the  calibration  range  in  order  to  correct  total 
temperature  measurements  (for  instance,  within  boundary  layers). 
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An  examination  of  the  methods  to  correct  total  temperature  mea- 
surements proposed  by  Bontrager  (Ref.  1)  and  others  (Refs.  4 and  5) 
was  undertaken.  An  analysis  was  made  of  the  differences,  similari- 
ties, and  problem  areas  of  the  three  correction  methods.  A new 
approach  to  the  total  temperature  probe  correction  was  developed  in 
an  attempt  to  account  for  the  shortcomings  of  these  previous  methods. 

The  analysis  is  based  on  the  total  temperature  variation  in  a laminar 
developing  flow  within  a tube  whose  walls  are  at  the  adiabatic 
recovery  temperature  of  the  local  flow  field.  This  approach  results 
in  the  ability  to  theoretically  correct  probe  data  for  all  local 
flow-field  conditions  with  calibration  data  only  required  in  the  free 
stream.  Corrections  made  to  total  temperature  measurements  by  this 
new  method  are  compared  to  theoretical  predictions  and  reveal  a sig- 
nificant improvement  over  previous  approaches. 

2.0  DISCUSSION  OF  PRESENT  METHODS 

The  total  temperature  probe  currently  in  use  at  AEDC-VKF  evolved 
from  the  need  for  a probe  that  would  give  accurate  and  consistent  results 
and  at  the  same  time  be  both  rugged  and  relatively  easy  to  manufacture. 
Figure  1 is  a sketch  of  the  probe  showing  the  important  physical  dimen- 
sions. Recent  tests  have  employed  probes  of  varying  size  depending  mostly 
upon  test  requirements.  Typical  configurations  that  have  been  used  are 
presented  in  Table  1. 

The  parametric  study  conducted  by  Bontrager  (Ref.  1)  was  used  to 
define  the  best  compromise  of  probe  geometry  in  order  to  reduce  the 
cumulative  errors  attributable  to  radiation,  velocity,  and  conduction. 

The  radiation  error  results  principally  from  the  inability  of  the 
probe  shield  to  block  energy  radiated  to  the  "cold"  tunnel  walls. 

Bontrager  examined  the  case  when  the  tunnel  flow  was  neither  an 
emitter  nor  an  absorber  and  concluded  that  the  radiation  error  could 
be  neglected  when  compared  to  the  convection  losses.  Velocity  error 
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Thermocouple 
Junction 


• Dv  Diameter  (2  Holes  180  deg  Apart 
Rotated  90  deg  from  Aft  Holes) 


Tube  Containing  Insulated 
Thermocouple  Leads 


Note-.  All  Dimensions  in  Inches 

Figure  1.  Schematic  of  single-shielded,  total  temperature  probe. 


Table  1 Typical  Values  of  Geometric  Parameters 
for  Probe  Shown  in  Fig.  1 


Configuration 

A 

D 

E 

°v 

J 

L 

1 

0.010 

0.034 

0.038 

0.012 

0.427 

0.068 

2 

0.020 

0.056 

0.060 

0.020 

0.668 

0.112 

3 

0.020 

0.076 

0.080 

0.027 

0.728 

0.152 

4 

0.040 

0.144 

0.156 

0.051 

1.032 

0.288 
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resulted  from  the  inability  of  the  thermocouple  junction  to  recover 
all  of  the  kinetic  energy  of  the  flowing  gas  within  the  probe.  For 
mean  probe  Mach  numbers  of  0.3  or  less,  the  velocity  error  was  shown 
to  be  less  than  0.25  percent.  Heat  conduction  along  the  thermoelectric 
wires  was  considered  by  modeling  the  thermocouple  wire  as  an  extended 
surface.  For  relative  high  length-to-diameter  ratios  of  the  wire, 
conduction  effects  were  shown  to  be  negligible. 

Bontrager's  analysis  indicated  that  probe  entrance  length-to- 
diameter  ratios  of  2.0,  entrance  area-to-vent  area  ratios  of  2 or 
greater,  and  exposed  wire  lengths  much  greater  than  the  wire  diameter 
would  result  in  total  temperature  measurements  essentially  free  of 
conduction,  radiation,  and  velocity  errors.  Probes  considered  here 
are  assumed  to  adhere  to  these  restrictions,  and,  thus,  the  effect 
of  these  errors  on  the  total  temperature  probe  correction  will  not  be 
considered  further. 


The  largest  single  factor  affecting  the  ratio  of  measured-to- 
actual  total  temperature  is  the  loss  in  energy  associated  with  the 
developing  thermal  boundary  layer  in  the  entrance  region  of  the  probe. 
Thus,  Bontrager  attemped  to  correlate  the  ratio  of  measured-to-actual 
temperature  with  the  products  of  Reynolds  number,  Re^,  based  upon 
conditions  at  the  probe  entrance,  where 


He,,  = 


7 

2 ^ 

\ * J_~ 

\ Poe” 

_v 

“TT 

) \r  - 1 

/ r!/,_ 

A'  PeO 

(1) 


and  the  ratio  D/L.  Here  A /A  is  the  vent  area-to-entrance  area 

v e 

ratio,  p0  is  the  pitot  pressure  at  the  total  temperature  probe  loca- 
6 

tion,  D is  the  probe  internal  diameter  (ID),  L is  the  distance  from 

the  probe  entrance  to  the  thermocouple  junction,  is  the  absolute 

viscosity  based  on  entrance  conditions,  and  Tn  is  the  local  total 

i 

temperature.  The  specific  heat  ratio,  y,  is  approximately  constant 
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: 


f 

i 


over  the  temperature  range  of  interest,  and  R is  the  gas  constant. 

The  functional  form  of  Re^  as  expressed  in  Eq.  (1)  assumes  choked  flow 

through  the  probe  vent  areas  (see  Fig.  1).  For  constant  Mach  number 

conditions,  calibration  data  appeared  to  correlate  Re^  with  the 

measured-to-actual  total  temperature  ratio  (T  / T ). 

o o . 
m 1 

The  main  error  in  Bontrager's  correlation  is  the  neglect  of  the 
loss  of  thermal  energy  associated  with  the  local  undisturbed  flow 
at  Mach  number  M^.  In  an  attempt  to  account  for  the  loss  in  kinetic 
energy  of  the  flow  when  brought  to  rest,  a recovery  factor  of  the 
probe  has  been  defined  and  used  in  numerous  studies  (for  example,  see 
Refs.  2 and  3).  For  a flow  brought  to  rest  adiabatically , 

To.  = Tj  (l  , ^ \lf)  (2) 

For  a nonadiabatic  flow,  a recovery  factor,  n,  may  be  defined,  in 
conventional  terms,  by 


where  T is  the  measured  total  temperature.  Eliminating  T.  from 
m 1 

Eqs . (2)  and  (3)  and  solving  for  n,  one  derives  Eq.  (A): 


Equation  (4)  has  been  used  extensively  as  a correlation  parameter 
against  Re^  for  a fixed  probe  geometry. 

Recently,  Laderman  (Ref.  4)  used  the  definition  of  n given  in 
Eq.  (4)  to  correlate  the  total  temperature  correction  with  the  local 
Reynolds  number,  Re^,  based  on  the  probe  diameter,  local  unit  mass 
flow,  and  local  total  temperature.  The  similarity  of  the  present 
shielded  probe  with  an  unshielded  thermocouple  probe  suggested  this 
particular  correction  procedure.  Even  though  this  method  does  attempt 
to  account  for  the  Mach  number  effect,  correlation  of  q with  the 
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local  Reynolds  number  overlooks  the  area  change  of  the  stream  tube 
captured  by  the  probe  because  of  changes  in  the  local  flow  condition 
(for  instance,  through  a boundary  layer).  Some  typical  calculations 
have  been  made  that  indicate,  for  a fixed  probe  geometry,  variations 
of  up  to  100  percent  in  the  diameter  of  the  stream  tube  captured  exist 
for  a probe  traverse  through  a 4-deg  cone  boundary  layer  at  Mach  10. 

By  a somewhat  similar  analysis,  Martellucci  (Ref.  5)  has  used  the 
unit  mass  flow  through  the  probe  versus  p for  correcting  total  tempera- 
ture measurements.  By  using  this  unit  mass  flow,  the  problem  encoun- 
tered in  Laderman's  method  is  eliminated.  The  use  of  unit  mass  flow  as 
a correlation  parameter  neglects,  however,  viscous  effects  in  the 
entrance  region  of  the  probe  for  length-to-diameter  (L/D)  values  dif- 
ferent from  zero.  Moreover,  for  a typical  flow-field  survey  as  men- 
tioned in  the  previous  paragraph,  the  absolute  viscosity  of  the  flow 
in  the  entrance  region  of  a probe  can  change  by  as  much  as  30  percent. 
Thus,  the  methods  of  Laderman  and  Martellucci  (Refs.  4 and  5,  respec- 
tively) effectively -deal  with  portions  of  the  correction  problem  but 
fail  to  handle  the  entire  problem  adequately. 

An  additional  serious  deficiency  that  pervades  all  the  previously 
mentioned  methods  is  the  fact  that  probe  calibrations  are  made  in  flow 
environments  considerably  different,  as  a general  rule,  from  the  envi- 
ronment over  which  the  corrections  are  applied.  Typical  calibrations 
of  total  temperature  probes  in  Tunnel  C generally  cover  an  Re^  range 
of  100  to  1,000.  Over  this  Re^  range  p is  well  defined.  For  values 
of  Re^  less  than  the  lower  limit  of  the  calibration  data,  considerable 
uncertainty  exists  in  the  shape  and  magnitude  of  the  n versus  Re^  curve. 
Moreover,  it  is  not  clear  that  the  definition  of  p properly  accounts 
for  the  Mach  number  effect. 

In  order  to  assess  the  magnitude  of  these  uncertainties,  some 
sample  calculations  were  made.  Calibration  curves  based  on  each  of 
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the  three  methods  were  allowed  to  vary  between  reasonable  limits  out- 
side the  range  of  the  calibration  data.  These  curves  were  then  applied 
to  the  correction  of  total  temperature  measurements  in  the  laminar 
boundary  layer  of  a 4-deg  sharp  cone  at  Mach  10  in  Tunnel  C.  Possible 
errors  in  the  corrected  measurements  of  10  to  100  percent  resulted. 
These  variations  are,  of  course,  extreme  and  nay  be  reduced  by  a more 
precise  definition  of  the  calibration  curve.  The  variations  are, 
however,  reasonable  in  light  of  present  calibration  procedures. 

3.0  IMPROVED  CORRECTION  METHOD 

The  previous  sections  have  examined  the  problems  associated  with 
three  methods  [those  of  Bontrager,  (Ref.  1),  Martellucci  (Ref.  5),  and 
Laderman  (Ref.  4)]  in  defining  an  accurate  correction  to  total  tempera- 
ture measurements  made  with  the  single-shielded,  thermocouple  probe. 
Deficiencies  have  been  shown  to  exist  in  these  methods  which  may  result 
in  large  errors  in  the  corrected  Tq  measurement.  The  purpose  of  this 
section  will  be  the  development  of  an  improved  correction  method  that 
effectively  eliminates  these  major  sources  of  error. 

Consider  the  flow  and  geometric  model  as  shown  in  Fig.  2.  Assume 

that  the  vent  flow  is  choked  over  the  entire  Mach  number  calibration 

and  correction  range.  The  developing  flow  in  the  entrance  region  of 

the  probe  is  assumed  to  be  laminar  for  all  entry  Reynolds  numbers 

based  on  the  distance  to  the  thermocouple  junction  and  to  be  unaffected 

by  thermal  gradients.  A developing  thermal  layer  exists  because  of 

the  heat  flux  resulting  from  the  difference  in  the  temperature  of  the 

entering  flow  (e)  and  in  the  adiabatic  recovery  temperature  of  the 

shield  based  on  local  conditions  (i).  The  size  of  the  thermocouple 

2 

junction  relative  to  the  internal  area  of  the  shield  (ttD  /4)  is 
assumed  to  be  small  enough  so  that  it  does  not  affect  the  flow  field 
within  the  tube.  Then,  by  defining  the  function  g as 
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where  T is  the  local  centerline  temperature  within  the  probe  and  where 

T and  T are  the  entering  temperature  of  the  flow  and  the  adiabatic 
e s 

recovery  temperature  of  the  shield,  respectively,  the  functional  form 
of  the  measured-to-actual  total  temperature  may  be  determined.  Hence, 


Here  c is  tue  specific  heat  at  constant  pressure.  The  terms  M and 
P e 

u/u  are  the  entrance  Mach  number  and  ratio  of  centerline  axial 
e 

velocity  at  the  thermocouple  junction  to  the  axial  entrance  velocity, 

respectively.  For  choked  vent  flow,  the  ratio  of  entrance  flow  area, 

A . to  effective  vent  flow  area,  A , defines  the  entrance  Mach  number 
e’  v 

for  isentropic  flow.  Equation  (6),  of  course,  implies  no  losses  in 
bringing  the  flow  to  rest  at  the  thermocouple  junction.  This  error 
is  negligible  for  all  practical  flow  geometries  considered. 


Figure  2.  Typical  probe  showing  parameters  of  importance. 
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Employing  the  conventional  definition  of  the  recovery  temperature  of 
the  shield,  one  may  refine  Eq.  (6)  to  give 


Tu 

1 + (,>  _ 1)  r<M;;  Mp,  h*) 


2cP\  e[w 


(7) 


where 


HMj;  M„.  R* ) 


I f 'Tj 


(1  + t- 


R*  M ; ) 


T.  1 

1 O 


(8) 


The  use  of  the  isentropic  flow  relations  allows  for  the  definition  of 

T /T . and  T . /T  . assuming  M and  M.  are  known.  The  recovery  factor 

of  the  shield,  R*,  is  equal  to  Prl/2  for  laminar  flow  where  Pr  is  the 

Prandtl  number  of  the  local  flow  field.  Figure  3 shows  the  influence 

of  the  recovery  factor  on  F as  a function  of  the  local  Mach  number. 

From  this  figure  it  may  be  seen  that,  for  fixed  M and  g,  increasing 

R*  tends  to  decrease  the  effect  of  the  local  Mach  number  on  the 

temperature  correction.  The  variation  of  T with  for  R*  equal  to 

0.84  is  shown  in  Fig.  4,  indicating  the  limiting  nature  of  r for 

increasing  values  of  A /A  . 

e v 


The  only  quantities  not  yet  defined  in  Eq.  (7)  are  the  function 
g and  the  velocity  ratio  u/ug.  It  is  appropriate  here  to  discuss 
the  possible  forms  of  g and  their  dependence  on  the  flow  field  within 
the  probe.  A comparison  of  various  solutions  will  be  made  first, 
followed  by  a brief  presentation  of  the  functional  form  of  the  solu- 
tions. Figure  5 shows  a total  of  five  solutions  for  g that  depend 
on  the  Peclet  number,  Pe,  which  is  equal  to  Re^Pr/2,  and  the  ratio 
L/D.  Generally,  for  all  values  of  L/D,  g is  bounded  by  the  values 
0 and  1,  asymptotically  approaches  these  limits  as  Pe  approaches  0, 
and  grows  without  bound.  A comparison  of  similar  curves  for  different 
L/D  values  indicates  a functional  shift  of  the  g curve  with  Pe.  An 
examination  of  the  effects  of  axial  heat  conduction  reveals  their  im- 
portance for  Peclet  numbers  less  than  20.  The  trends  of  g versus 
Peclet  number  for  comparable  fully  developed  laminar  and  potential 
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1 


(Recovery  Factor)  R* » 0. 80 


showing  el 
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flows  also  suggests  a functional  shift  of  the  g curve  for  a given  L/D 

value.  The  variation  of  g with  Pe  and  L/D  was  evaluated  numerically 

for  a developing  flow  in  the  entrance  length  of  the  probe  for  the 

case  with  no  axial  conduction.  In  the  limits,  the  developing  flow 

solution  tends  to  fully  developed  laminar  and  potential  solutions  for 

Pe  approaching  0 and  for  Pe  approaching  infinity,  respectively. 

Included  in  these  figures  is  a transition  curve  derived  by  shifting 

the  developing  flow  solution  an  amount  equal  to  the  difference  in  the 

potential  flow  solutions  with  and  without  axial  conduction.  Thus, 

the  required  g curve  to  be  used  in  Eq.  (7)  for  the  total  temperatuie 

correction  is  the  transition  curve  coupled  with  the  developing  flow 

solution  for  u/u  . 

e 

The  five  basic  solutions  were  derived  from  the  governing  energy 
equation  for  axisymmetric , steady  flow  in  the  axial  direction  only. 

The  effect  of  the  radial  velocity  component  in  the  entrance  region  of 
the  probe  on  the  thermal  layer  solution  is  neglected  since  Pe  is 
generally  less  than  300  (Ref.  6).  The  effect  of  the  axial  velocity 
gradient  in  the  radial  direction  has  also  been  assumed  negligible. 

The  boundary  conditions  and  governing  equation  in  terms  of  the  func- 
tion g and  a)  = r/rQ  and  x'  = x/(rQPe)  are  given  in  Eqs.  (9)  and  (10) 
with  coordinates  as  illustrated  in  Fig.  2;  here  r^  is  the  tube  radius 
equal  to  D/2. 

dg  J d2g  ]__$_(  §l\  (9) 

p,2  a*'2  V H 

e 

B.  C.:  g = 0 for  x ' > 0,  (o  - 1 

g = 1 for  x ' = 0,  0 < at  < 1 

g 0 as  x'  » (10) 

An  exact  solution  to  Eq.  (9)  subject  to  the  boundary  conditions  given 
in  Eq.  (10)  is  possible  for  the  case  where  u/ug  = 1.0  (potential  flow). 
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The  solution  may  be  expressed  in  the  form  of  Eq.  (11)  for  the  cases 
considered,  neglecting  axial  conduction  represented  by  the  term 


K 


2 

anJl<an> 


2 

exp  J0(anaj) 


(ID 


where 


neglecting  axial  conduction  (Ref.  7)  and 


(12) 


(13) 


including  axial  conduction  (Ref.  8).  and  are  Bessel  functions 

of  the  first  kind  of  order  zero  and  one,  respectively,  and  are  the 
zeros  of  J^a^).  An  approximate  analytic  solution  for  the  fully 
developed  laminar  case  including  axial  conduction  (Ref.  8)  is  also 
shown  to  be  expressed  by  Eq.  (11)  with,  however. 


ft 


2 

n 


(14) 


2 

where  u/u  = 2(1  - w ).  It  is  important  to  note  that  the  three  solu- 
e 

tions  can  all  be  expressed  in  the  form 


g = 


i C„ 

n=  1 " 


exp 


(15) 


for  w = 0,  where  the  C are  constants  and  q is  a function  of  Pe  and 

n n 

cx^  only.  This  similarity  reinforces  the  observation  made  earlier 
that  the  effect  of  a given  change  in  L/D  is  to  shift  the  g curve  with 
respect  to  Peclet  number  by  approximately  a constant. 


The  numerical  evaluation  of  the  case  for  developing  flow  neg- 
lecting axial  conduction  effects  follows  Kays  (Ref.  9),  employing, 
however,  the  solution  of  Sparrow,  Lin,  and  Lundgren  (Ref.  10)  for  the 
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hydrodynamic  flow  development  in  the  entrance  region  of  the  tube. 

The  axial  velocity  distribution  within  the  tube  is  given  in  Eq.  (16): 


2(1  - 


where  x*  is  evaluated  by 


2 - _4  J0(tVJ> 

' * n'l  J0<an) 

n L 


f’rx  ' = / ( d 


O 

exp  ( -a~  x*) 


and  e is  given  by 


>r  .-Kt) 

l *£)-  '■*£)  ^ 


r)(ij  Ifxj  - I Q L doj 

Equation  (18)  has  been  evaluated  numerically  and  is  given  in  Fig.  6 
for  x*  < 0.20.  For  values  of  x*  greater  than  0.20,  £ is  approximately 
constant  and  equal  to  1.82.  Figure  7 is  a plot  of  the  axial  velocity 
ratio  along  the  tube  centerline  as  a function  of  Re^/ (L/D) . 

Using  the  expression  for  the  velocity  ratio  given  in  Eq.  (16), 
the  governing  energy  equation  was  solved  numerically  for  the  case  of 
no  axial  conduction.  For  the  interior  grid  points 


*l.i 


where 


0(£)[i  ‘ i] 

•• 

■ mi ^ - i] 
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Here  Ax'  and  Aw  are  the  grid  spacings  in  the  axial  (i)  and  radial  (j) 
directions,  respectively.  Along  the  tube  centerline. 


Bi+l.l  ’ f2  I 

and  at  the  boundaries. 


(21) 


g = 0 at  gj  = 1,  x ' > 0 
p = 1 at  x - 0,  0 ^ (u 


(22) 


Typical  velocity  and  temperature  profiles  are  shovm  in  Fig.  8.  The 
general  solution  for  g as  a function  of  Pe  and  L/D  along  the  tube 
centerline  is  shown  in  Fig.  9 with  the  axial  conduction  effect  in- 
cluded. Thus,  with  the  definition  of  u/u^  and  g as  a function  of 
Pe,  ReQ,  and  L/D,  an  improved  calibration  and  correction  of  total 
temperature  probe  data  may  be  obtained  through  the  use  of  Eqs.  (7) 
and  (8). 


4.0  DISCUSSION  AND  CONCLUSIONS 

The  preceding  section  has  presented  a new  method  for  the  cor- 
rection of  total  temperature  measurements  made  with  a single-shielded, 
thermocouple  probe.  A discussion  of  some  of  the  important  character- 
istics and  their  relationship  to  previous  methods  (Refs.  1,  4,  and  5) 
will  aid  in  an  understanding  of  the  applicability  of  the  new  technique. 

An  application  of  the  new  method  to  some  laminar  flow  data  from 
a recent  test  in  Tunnel  C is  given  in  Figs.  10  through  13.  The  total 
temperature  probe  used  had  nominal  dimensions  corresponding  to  config- 
uration 2 as  listed  in  Table  1.  The  experimental  calibration  data 
were  obtained  in  the  free  stream  by  varying  the  tunnel  stagnation 
pressure.  The  data  were  used  to  determine  the  experimental  value  of 
the  vent  ratio.  Figure  10  shows  the  fit  of  the  g curve  through  the 

calibration  data  for  the  experimental  value  of  A /A  of  0.11.  The 

v e 

experimental  data  scatter  about  the  g curve  is  equivalent  to  a maximum 

error  in  T0  / TQ  of  the  theory  compared  to  the  calibration  data  of  ±0. 5 
m i 

percent. 
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Figure  9.  Temperature  distribution  along  probe  centerline  as 
a function  of  Peclet  number  for  developing  flow 


with  axial  conduction. 


I 


i.  Calibration  curve  for  improved  method. 
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Figure  11.  Comparison  of  the  corrected  experimental  and 
theoretical  total  temperature  within  the  bound- 
ary layer  of  a 4-deg  cone  at  low  Re  . . 
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Figure  13.  Comparison  of  the  corrected  experimental  and 
theoretical  total  temperature  within  the  bound 
ary  layer  of  a 4-deg  cone  at  high  Re^L . 
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In  general , the  vent  area  ratio  is  a function  of  the  pressure 

ratio,  p./p0  , and  the  vent  Reynolds  number.  Re  (see  Ref.  11).  This 
i 0 v 

effect  has  been  included  in  the  present  correction  procedure.  For  the 

calibration  and  flow-field  data  considered  here,  the  variation  of  A /A 

v e 

with  Re^  and  p^/p0  was  insignificant.  For  smaller  probes,  however, 
this  effect  will  become  important.  A post-test  bench  calibration  of 
the  probe  gave  A^/A^  = 0.14.  The  low  experimental  value  as  compared 
to  A^/A^  = 0.5  obtained  from  Table  1 is  caused  by  the  partial  restric- 
tion of  the  vent  holes  during  construction.  The  character  of  the  mean 
flow  over  the  vent  holes  was  expected  to  be  a significant  parameter  in 
determining  the  effective  vent  flow  area.  The  difference  in  the  value 

of  A /A  obtained  by  calibration  in  the  tunnel  and  by  the  bench  cali- 
v e 

bration  (correcting  for  relative  Re^  and  p^/pog)  was  felt  to  be 
representative  of  this  effect  and  the  difference  was  assumed  constant 
for  the  local  flow-field  conditions  considered. 

Measurements  were  taken  at  an  axial  station  of  0.88  of  the  model 
length  measured  from  the  nose  on  the  windward  ray  of  a 4-deg  sharp 
cone  at  a nominal  angle  of  attack  of  1 deg  and  Mach  number  of  10.  All 
distances  referred  to  in  Figs.  11,  12,  and  13  have  been  nondimension- 
alized  by  the  model  axial  length. 

Figure  1 1 shows  the  uncorrected  and  corrected  total  temperature 
data  as  a function  of  the  Z'  coordinate  measured  normal  to  the  cone 
surface  at  a Reynolds  number  based  on  total  model  length  (Re^  ) of 

0.588  x 10^  and  wall  temperature  ratio  of  T /T  = 9.96.  The  band  on 

w 00 

the  corrected  data  represents  a probe  location  uncertainty  of  ±0.0006 

in  Z'.  For  comparison  purposes,  the  hypersonic  viscous  shock  layer  (HVSL) 

theory  of  Ref.  12  is  also  shown  in  Fig.  11  for  identical  free-stream 

conditions.  The  theoretical  total  temperature  profiles  were  computed 

using  an  "effective"  c based  on  the  average  of  the  local  static  and 

P 

stagnation  temperatures.  The  corrected  data  are  seen  to  be  in  excellent 
agreement  with  the  theory.  This  is  indirect  verification  of  the  present 
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I 


method  since  there  were  large  variations  within  the  boundary  layer  of 
Re^  (from  1.5  to  50)  and  M.  (from  1.2  to  9.82)  compared  to  the  calibra- 
tion conditions  of  23  < Re  < 120  and  M = 10.  Also  indicated  in  the 

D ® 

figure  is  the  relative  size  of  the  probe.  Considering  the  probe  size 
effect  near  the  wall  and  at  the  knee  in  the  T^/T^  profile,  the  corrected 
data  also  exhibit  the  proper  trends. 

Figures  12  and  13  give  a comparison  of  the  corrected  data  and 
theory  for  Re  , = 1.85  x 106,  T /T  = 10.74  and  Re  , = 2.74  x 106, 

T /T^  = 11.72,  respectively.  Agreement  here  is  seen  to  be  good  in  both 
magnitude  and  expected  trend.  The  effect  of  probe  size  on  the  measured 
data  for  these  higher  Reynods  number  conditions  when  compared  to  Fig.  11 
is  greater  at  the  knee  of  the  total  temperature  profile  because  of  the 
steeper  temperature  gradient.  Thus,  the  present  correction  results  in 
very  good  total  temperature  correlation  between  experiment  and  theory. 


An  interesting  aspect  of  the  new  correction  method  is  its  relation- 
ship to  the  methods  of  Laderman  and  Martellucci  (Refs.  4 and  5,  respec- 
tively). If  the  total  temperature  correction  equation,  Eq.  (7),  is 
rewritten  using  Eq.  (4)  with  n as  defined  in  Refs.  4 and  5,  the  correction 
equation  may  be  reduced  to  the  following  form: 


V 


(g-  DP'  + 


(23) 


where 

P'  = r/(l  - T/TJ  (24) 

1 

For  very  large  values  of  Ag/Av  (on  the  order  of  10  or  greater)  and 
local  Mach  numbers  much  larger  than  1,  the  contribution  of  the  term 
representing  kinetic  energy  of  the  developing  flow  within  the  tube 
in  Eq.  (23)  may  be  neglected  and  T'  is  approximately  constant.  Thus, 
n becomes  a function  of  ReD  and  L/D  only.  This  result  agrees  with  the 


.i. 
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form  of  n used  by  Laderman  and  Martellucci.  As  the  entrance  area-to- 
vent  area  ratio  is  decreased,  however,  the  kinetic  energy  term  increases 
in  magnitude  and  may  not  be  neglected.  Moreover,  ?'  may  no  longer  be 
assumed  independent  of  M^. 

In  some  of  Bontrager's  experimental  results,  there  is  a definite 

''overshoot"  of  the  calibration  curve  above  1.0.  Bontrager  could  not 

account  for  the  magnitude  of  this  overshoot.  The  form  of  the  new 

correction  as  given  in  Eq.  (7)  does  allow  for  an  overshoot  in  TQ  /T0 

m i 

that  is  predominantly  influenced  by  the  value  of  A^/A^.  For  example, 

under  free-stream  calibration  conditions  (Re^  > 200)  with  L/D  = 2, 

values  of  A^/A^  on  the  order  of  two  are  required  to  give  a two-percent 

overshoot  in  TQ  /TQ  . Another  feature  of  the  new  method  when  compared 
m i 

to  the  previous  methods  applied  to  the  calibration  and  correction  of 

total  temperature  probe  data  is  the  fact  that  the  new  method  has  a 

lower  limit  on  TQ  /Tq  . For  a given  probe  configuration,  the  limit 
m i 

is  approximately  equal  to  1 - T for  large  values  of  Ag/Av  and  occurs 

when  T = T . 

e 

In  conclusion,  methods  of  total  temperature  correction  for  a 
single-shielded,  thermocouple  probe  have  been  examined.  The  important 
problems  associated  with  each  method  have  been  addressed.  Because  of 
the  inability  of  any  one  of  the  earlier  methods  to  handle  satisfac- 
torily the  correction  problem,  a new  method  was  developed  based  on  a 
detailed  examination  of  the  developing  thermal  layer  within  the  tube. 

The  theoretical  form  of  the  solution  allows  for  total  temperature 
corrections  to  be  made  in  an  uncalibrated  Mach  number  and/or  local 
Reynolds  number  region.  This  new  correction  method  appears  to  be  a 
significant  improvement  over  previous  approaches. 
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NOMENCLATURE 

A Flow  area 

C Series  constants 

0^  Specific  heat  at  constant  pressure 

D Probe  entrance  diameter 


31 


1 


w 


AEDC  TR  76  140 


D 

V 


f rf2,f3 


M 


n 


Pe 


Pr 


q 

R 


R* 

Re 

Re 


D 

D 


Re 

v 


Probe  vent  hole  diameter 
Functions  as  given  in  Eq.  (20) 

Function  equal  to  (T  - Tg)/(Te  - Tg) 

Entrance  length  to  thermocouple  junction 

Mach  number 

Integer 

Peclet  number  equal  to  PrRe^/2 

Prandtl  number 

Pressure 

Total  pressure 

Function  of  Pe  and  a 

n 

Gas  constant 

Recovery  factor  of  the  shield 

Reynolds  number  based  on  probe  entrance  condition 

Reynolds  number  based  on  local  condition  and  probe 
entrance  diameter 

Probe  vent  Reynolds  number 
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Re 


,L 


r 


r 

o 


T 


T 

o 


u 


X 


x' 


X* 


Z ' 


a 


a 

n 


r 

r' 

Y 


Free-stream  Reynolds  number  based  on  cone  model  length 

Radial  coordinate  as  shown  in  Fig.  2 

Probe  entrance  radius 

Temperature 

Total  temperature 

Axial  velocity 

Axial  coordinate  as  shown  in  Fig.  2 
Dimensionless  axial  coordinate 
Stretched  dimensionless  axial  coordinate 
Dimensionless  surface  normal  coordinate 
Angle  of  attack 

Eigenvalues  equal  to  zero's  of  J^a^) 

Function  of  a and  Pe  as  given  in  Eqs.  (12),  (13), 
and  (14) 

Function  as  defined  in  Eq.  (8) 

Function  as  defined  in  Eq.  (24) 

Specific  heat  ratio 
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Ax  ' 


Auj 


e 


n 


u 


u) 


SUBSCRIPTS 


e 


i 


s 


v 


w 


m 


Grid  spacing  in  x'  direction 

Grid  spacing  in  oj  direction 

Stretch  function  as  given  in  Eq.  (18) 

Total  temperature  recovery  function  as  given  in  Eq.  (A) 
Absolute  viscosity 

Dimensionless  radial  coordinate  equal  to  r/rQ 


Conditions  at  probe  entrance 


Local  undisturbed  conditions 


Conditions  at  shield  surface 


Conditions  at  shield  vent 


Conditions  on  the  cone  surface 


Conditions  in  the  free  stream 


Measured  conditions 
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